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Abstract 

We show that the hypothesis of abelian dominance allows QCD-monopoles 
to preserve the topological feature of the QCD vacuum within SU(2) lattice 
gauge theory. An analytical study is made to find the relationship between 
the topological charge and QCD-monopoles in the lattice formulation. The 
topological charge is explicitly represented in terms of the monopole current 
and the abelian component of gauge fields in the abelian dominated system. 
We numerically examine the relation and demonstrate the abelian dominance 
in the topological structure by using Monte Carlo simulation. 
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I. INTRODUCTION 



It is important to understand topological aspects of QCD to interpret several basic prop- 
erties of the vacuum structure. Two distinct pictures of the QCD vacuum, which are based 
on the presence of topological objects; instanton and QCD-monopole, are now widely ac- 
cepted. The main reason is that they give nice descriptions of several non-perturbative 
features, e.g. chiral symmetry breaking and color confinement. As is well known, SU(iV c ) 
gauge theory has classical and non-trivial gauge configurations (instantons) satisfying the 
field equation in the 4-dimensional Euclidean space R 4 PJ. The topological charge, which 
corresponds to the non-trivial homotopy group 7r3(SU(iV c )) = Z^, is assigned to instanton 
configurations jy. It has been established that such a topological feature plays an essential 
role on the resolution of the Ua(1) problem ||. Furthermore, the instanton liquid char- 
acterized by a random ensemble of instanton and anti-instanton configurations succeeds in 
explaining several properties of light hadrons, e.g. spontaneous chiral-symmetry breaking 
(S X SB) 0. 

Another picture of the QCD vacuum is motivated from the stimulating idea of abelian 
gauge fixing, which was proposed by 't Hooft and also independently by Ezawa-Iwazaki 
After performing a partial gauge fixing, which leaves an abelian gauge degree of free- 
dom, one knows that point-like singularities in the three-dimensional space R 3 , associated 
with the homotopy group 7r 2 (SlJ(A^ c )/U(l) 7Vc ~ 1 ) = -Z^ c_1 , can be identified as magnetic 
monopoles (QCD-monopoles) If monopoles are condensed in the QCD vacuum, the 

dual Meissner effect results f|. Recent lattice QCD simulations actually show that the dual 
Meissner effect is brought to the true vacuum by QCD-monopole condensation^ (see, e.g. a 
recent review article ||). Hence, color confinement can be regarded as the dual version of 

^QCD-monopole condensation is characterized by the presence of the long and tangled monopole 
trajectories in the four-dimensional space R 4 and can be interpreted as the Kosterlitz-Thouless 
type phase transition. 
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superconductivity. 

It seems that instantons and QCD-monopoles are relevant topological objects for the 
description of distinct phenomena. Here, we should emphasize that QCD-monopoles would 
play an essential role on non-perturbative features of QCD, which include not only confine- 
ment and but also S^SB. This possibility was first studied by using the Schwinger-Dyson 
equation with the gluon propagator in the background of condensed monopoles |7j . The idea 
of providing S^SB due to QCD-monopole condensation was supported by lattice simulations 
fHH. Thus, these results shed new light on the non-trivial relation between instantons and 



QCD-monopoles. Recently, both analytic works |[10| , pT| and numerical works |I^Jl^JlqJi4| . pl3 
have shown the existence of the strong correlation between these topological objects in spite 
of the fact that they originate from different homotopy groups. Furthermore, monopole tra- 
jectories become more complicated with the instanton density increasing in the background 



of a random ensemble of instanton solutions ]T7|. These results seems to suggest that the 



instanton liquid and QCD-monopole condensation may be indistinguishable. 

In this paper, we would like to know whether the correlation between instantons and 
QCD-monopoles has some physical significance, or not, from the topological viewpoint. 
Here, we must not forget the hypothesis of abelian dominance, which Ezawa and Iwazaki 
had first advocated 0. This hypothesis is essentially composed of two sentences Q: 

• Only the abelian component of gauge fields is relevant at a long-distance scale. 

• The non-abelian effects are mostly inherited by magnetic monopoles. 

Actually, lattice Monte Carlo simulation indicates that abelian dominance for some physical 
quantities, e.g. the string tension |r8|JT9| , the chiral condensate and also several light 



hadron spectra pO| , |2T| is realized in the maximally abelian (MA) gauge as well as monopole 
dominance. In this gauge, at least, the abelian component of the gauge field could be an im- 
portant dynamical degree of freedom at a long-distance scale. Here, an unavoidable question 
arises relating to the non-trivial correlation between instantons and QCD-monopoles. In the 

abelian dominated system, is it possible for the non-abelian topological nature to survive? For 



such an essential question, Ezawa and Iwazaki have proposed a remarkable conjecture p2| : 
once abelian dominance is postulated, the topological feature is preserved by the presence of 
monopoles. The main purpose of this paper is to confirm the presence of the non-trivial cor- 
relation between instantons and QCD-monopoles by finding evidence for the Ezawa-Iwazaki 
conjecture. 

For simplicity, we restrict the discussion to the case of SU(2) gauge group throughout this 
paper. The organization of our paper is as follows. In Sec. 2, we show that the topological 
charge is explicitly represented in terms of the monopole current and the abelian component 
of gauge fields through the hypothesis of abelian dominance. In Sec. 3, we confirm numeri- 
cally its justification by means of Monte Carlo simulation within SU(2) gauge theory after 
the MA gauge fixing. Finally, Sec. 4 is devoted to a summary and conclusions. 

II. ABELIAN DOMINANCE HYPOTHESIS FOR TOPOLOGICAL CHARGE 



We first address the definition of the topological charge in the lattice formulation. The 
naive and field-theoretical definition of the topological density |23j is given by 



q{s) = ^e fjLUpa ti{P^(s)P pa {s)} , (1) 

with the clover averaged SU(2) plaquette: 

P^s) = i ( U^(s) + Ul^s) + Ul_ v {s) + U^_ v {s) ) . (2) 

Here we have used the convenient notation for the SU(2) link variable; £/_^(s) = UUs—fi) so 
that U± p ± u {s) = U± Pj (s)U± u (s+jl)U± fJj (s+C')U± pj (s). One naively expects that the topological 
charge Q CO nt is extracted from the summation of the previously defined topological density 



over all sites up to leading order in powers of the lattice spacing a [23 



Ql = £ - Qcont + O(a') , (3) 

where Qcont = S s tr {a 4 g 2 G^(s)*G M!/ (s)}. Strictly speaking, the value Ql has not only 
0(a 6 ) corrections, but also a renormalized multiplicative correction of Q con t i n Eq-© ||24j| . In 



general, one needs some smoothing method to eliminate undesirable ultraviolet fluctuations 
in the numerical simulation so that a renormalization factor would be brought to unity. 

If we assume that the QCD vacuum is described as the abelian dominated system in a 
suitable abelian gauge, then the SU(2) link variable is expected to be U(l)-like as U^(s) ~ 
u p (s) = exp {ia 3 9^(s)}. Here, we define the angular variable 9^ as 

^( S )=arctan[^ 3 ( S )/^( S )] (4) 

in the compact domain [— ir, ti) when the SU(2) link variable is parameterized as = 
^°( s ) + ^ a aU^(s) (a = 1,2,3) pjfl . In the abelian dominated system, we might consider 
the abelian analog of the topological density q Abel (s) flEfl , instead of the previously defined 
topological density, through the replacement of P pv by the clover averaged U(l) plaquette 

Pflv- 

P^u(s) = - ( ?v(s) + u ] _ pv (s) + ul_ u (s) + u^^ u (s) ) 

1 1 

= tE u^s - ijl - ji>) , (5) 

4 i,j=0 

where u py denotes the U(l) elementary plaquette. The explicit expression of q Abcl (s) is then 
given by 

9 A bci( s ) = 2 e M^ tr {P^i.s)p pu (s)} 

1 1 

= -— £ fiu P aSin9 fiu (s -ifi- jD)sm9 pa {s -kp- la) , (6) 

ib i,j,k,l=0 

where 9^ v {s) = d p 9 u (s) — d u 9 p (s) p6| . <9 M denotes the nearest-neighbor forward difference 
operator satisfying d^fi^s) = f(s + p) — f(s). 

Our next aim is to discuss the expression of the abelian analog of the topological density 
in the naive continuum limit a — > 0. This is because we need only the leading order term in 
powers of the lattice spacing to determine the corresponding topological charge. Here, one 
may notice that the U(l) elementary plaquette is a multiple valued function of the U(l) 
plaquette angle 9^ v . Hence, we divide 9^ into two parts, 
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M s ) = W s ) + 27rn^(s) , (7) 

where is defined in the principal domain [—71,11), which corresponds to the U(l) field 
strength. The anti-symmetric tensor Thiols CcLll take the restricted integer values 0, ±1,±2. 
Taking the limit a — > 0, i.e. Q pv — > 0, we thus arrive at the following expression [2B|: 



?Abd( s ) ~ -W e ^(s)9p^(«) , (8) 

where 0^(s) = \ J2lj=o — if 1 — jv)- The r.h.s of Eq.® does not necessarily vanish in 
spite of the fact that it is represented only in terms of the U(l) field strength. 

Next, we show that the non- vanishing contribution to the value of <Z Abel (s) results from 
monopoles. For the identification of monopoles, we follow DeGrand-Toussaint's definition 
in the compact U(l) gauge theory ||27|| . The magnetic current is given by 

1 

k^{s) = —e pupa d u 9 pa (s + ft) 

47T 

= -]^^ P ud u n pa {n + p) . (9) 

In the last line, we have used the Bianchi identity on the U(l) plaquette angle; e pi , pa d u 6 pa = 0. 
Then we can easily see that the current k^ denotes the integer- valued magnetic current, i.e. 
the monopole current [27]. Strictly speaking, the monopole current resides at the dual lattice 



site orthogonal to direction fi. For simplicity, we have used an undiscriminating notation 
between the ordinary lattice site and the dual lattice site. One may notice that it is obvious 
that k^ is topologically conserved; d'^k^ = 0. & denotes the nearest-neighbor backward 
difference operator satisfying & f(s) = f(s) — f(s — ft). 

Here, we define the averaged magnetic currentj] in the 2 3 extended cube at the center of 

as 



the lattice site s orthogonal to direction /i [p£ 

1 



M s ) = g £ k p (s-iv - jp-ka) , (10) 

i ,j,k=0 



Tin some sense, this is similar to the type-II extended monopole current (2 
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where there indices (v,p,o~) are complementary to p. Using the nearest-neighbor central 
difference operator; A /t = \{d p + & }, we therefore rewrite as the following form: 

^m( s ) = ^vv P aA u Q pa (s) , (11) 

which obviously satisfies the conservation law; A M /C M (s) = 0. 

To show the explicit contribution of monopoles to the topological charge, we introduce 
the dual potential £> M satisfying the following equation p| : 



(a 2 ,^ - A M A„) = -2tt/C m ( S ) . (12) 

where A 2 = A M A M . We can perform the Hodge decomposition on the clover-average U(l) 
field strength with the dual potential B p as 

e^(s) = A M A(s) - A^s) + e^ pa A p B a {s) , (13) 



where A^ is the Gaussian fluctuation p9| . After a little algebra using the partial summation, 
we find the explicit contribution of monopoles to the r.h.s of Eq.(p|) |26[ as 



&pvpo®pv (s) @p<r (s) 4£^ / \ u ,£ - \ t<; p ( j A^vAy (s) Api3(j(s) 

= 16vr^( S )/C M ( S ) + • • • , (14) 

where the ellipsis stands for the total divergence, which will drop in the summation over all 
sites. Consequently, we arrive at the conjecture that the topological feature is preserved by 
the presence of monopoles in the abelian dominated system |[26|| : 

Qcont ^ - J2 V„ ( s ) , ( 15 ) 

where q Mono (s) = — 167r^4 A1 (s)/C A1 (s). In addition, the Gaussian fluctuation can not be defi- 
nitely identified except for the Landau gauge solution; A^A^s) = where a superscript L 
denotes 'Landau'. However, we need only to know the Gaussian fluctuation in the Landau 
gauge on the measurement of the quantity J2 S 9moho( s )> which has the residual U(l) gauge 
invariance. 



III. NUMERICAL RESULTS 



This section is devoted to a numerical analysis to justify the above conjecture through 
Monte Carlo simulation within SU(2) lattice gauge theory using the standard Wilson action. 
We choose the MA gauge as an applicable abelian gauge for the realization of abelian 
dominance in this paper. As we pointed out, the QCD vacuum has been observed as an 
abelian dominated system by lattice Monte Carlo simulations ||. We shall return to the 
explicit procedure of the MA gauge fixing later. 

As we have mentioned before, one needs to smooth the Monte Carlo gauge configurations 
in order to determine the topological charge. To eliminate undesirable fluctuations, we adopt 
the naive cooling method which is an iterative scheme to replace each link variable using 



the following procedure |23[ : 



U,(s)^U'^s) = ^(s)/\\J:,(s)\\ , (16) 
where HE^s)!! denotes the square root of determinant of £ M and 

= £ {U v {s)U,{s + v)Ul{s + fi) + Ul(s - v)U,{s - u)U v {s + fi - v)) . (17) 
This procedure is derived from the condition to minimize the action through the variation 



of the link variable ||23|| . Then, an iterated process leads to a local minimum of the action, 
which corresponds to a solution to the field equation. The relation Q L « Qcont is expected 
through the cooling procedure. 

Before turning to the abelian gauge fixing, we must take account of the permutation 
between the cooling procedure and the gauge fixing procedure. The action in the compact 
lattice gauge theory is not usually affected by the gauge fixing procedure for the link variables 
because of the gauge-invariant measure. However, the fundamental modular region |30 



appears in the case of the MA gauge fixing ||31|| . Then, the action is modified by the 
contribution from the Faddeev- Popov determinant [[3TJ. In other words, the above cooling 
procedure is not justified after the MA gauge fixing. Thus, we have to apply the MA gauge 
fixing procedure to the resulting gauge configuration after several cooling sweeps. 
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In the lattice formulation, the MA gauge fixing is defined by maximizing the gauge 
dependent variable R under gauge transformations; U^(s) — ► U^(s) = £l(s)Un(s)£l>(s + jl) 



R[n}=Y,tT{a 3 U,(s)a 3 Ul(s)}=2j2(l-2 X ,(s)) , (18) 

where %m = (^( s )) 2 + {U^{s)) 2 . The maximization of R implies that the off-diagonal 
components of the SU(2) link variables are minimized as much as possible through the 
gauge transformation. 

After the MA gauge fixing, we extract the U(l) field strength 9^ v and the monopole cur- 
rent from the U(l) link variable following the previous section. The Gaussian fluctuation 
in the Landau gauge is computed by the convolution of the electric current with the lattice 
Coulomb propagator T>(s — s') 

A^s) = -J2V(s-s')A x e x ,(s>), (19) 

s' 

where the lattice Coulomb propagator satisfies the equation; A 2 T>(s — s') = —5 S , s >. One can 
obtain the explicit form of T>(s — s') on a L 4 lattice as 

Vis - s') = V - , e ip < s ~ s,) (20) 

1 ; VE^sin 2 ^) 1 ; 

with the following abbreviation: = ^j-n^ and J2 P = jr^-^Hn~=o- 

We generate the gauge configurations by using the heat bath algorithm on a 16 4 lattice at 
(3 = 2.4. All measurements have been performed on independent configurations, which are 
separated from each other by 500 sweeps after the system has reached thermal equilibrium 
for 1000 heat bath sweeps. For the MA gauge fixing in the numerical simulation, we use the 
overrelaxation method with parameter u — 1.7 [0. Let us concentrate on the realization 
of the abelian dominance for topological density in the MA gauge. Before turning to the 
topological density, we show the probability distribution of the amplitude of the off-diagonal 
component Xn m Fig.[TJ. The solid curve, the broken curve and the dotted curve denote 
the case for configurations before cooling and after cooling for 1 sweep and 30 sweeps, 
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respectively. This figure tells us that the off diagonal components of the SU(2) link variable 
are actually made as small as possible on the whole lattice in the MA gauge. Furthermore, 
such a tendency becomes prominent with more cooling sweeps. The information from Fig.[L] 
suggests that large numbers of cooling sweeps could enhance abelian dominance. In addition, 
the off diagonal components; {7* and behave like random variables on the whole lattice 
before the MA gauge fixing so that the probability distribution is flat on unity in the case 
of no gauge fixing whether configurations are cooled, or not. 

We have checked the above expectation on the measurement of the topological density 
q(s) and the abelian analog of topological density gAbei(s). Fig. 2 and Fig. 3 display q(s) 
and (?Abei( s ) in a two-dimensional plane after 30 and 100 cooling sweeps. We define this 
plane as a slice at the center of the instanton in the (z, i)-plane. The center is identified as 
the local maximum of q(s). The topological density around the instanton is little affected 
by the cooling procedure as is shown in Fig. [2(a)| and Fig. |2(b)| . For the abelian analog of 



topological density, such a stability in the cooling process is not required because of the lack 
of the topological basis in the U(l) manifold. At 30 cooling sweeps, the abelian dominance 
for topological density is not largely revealed as shown in Fig. |3(a)| , since a large amplitude 
of the off-diagonal components still remains locally, especially around the center of the 
instanton and/or the monopole. However, after 100 cooling sweeps, one can recognize that 
a similar lump to the topological density is located around the center of the instanton in 
Fig. |3(b)| . Thus, the topological density is dominated by the abelian analog of topological 
density after enough cooling sweeps, as expected. 

Next, we will examine the correlation between topological charge and monopoles. The 
two types of corresponding topological charge are computed through the following formula: 

Qsuw^-^EsOO. ( 21 ) 

Q Mo no = -Jq^ E Vo ( s ) • ( 22 ) 

Qsv(2) denotes the ordinary topological charge. Qmoiio denotes the corresponding topological 
charge via monopoles. We show the scatter plots of Qsu(2) vs. QMono in Fig.4; (a) no 
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gauge fixing and (b) the MA gauge fixing after 100 cooling sweeps for 50 independent 
configurations. Obviously, there is not any correlation between the two topological charges 
before the MA gauge fixing. A one-to-one correspondence between Qsu(2) and Qmoiio is 
revealed in the scatter plot after the MA gauge fixing. In further detail, we find that there 
is a small variance between Qsv(2) and Qmoiio- The slope of the linear correlation in the 
scatter plot is not unity, but 1.43 in Fig. [4(b)| (16 4 at (3 = 2.4). On the several measurements 
at (3 = 2.45 and 2.5, it seems that this slope hardly depends on the lattice spacing after 
large enough numbers of cooling sweeps (see Table 1). Table 1 tells us that the relation; 
QMono ~ 0.7Qsu(2) is almost satisfied in the MA gauge on these data. 



p 


cooling sweeps 


30 


50 


100 


2.40 


1.45 


1.46 


1.43 


2.45 


1.46 


1.43 


1.44 


2.50 


1.49 


1.43 


1.43 



To discuss the topological feature on QMono, we show the probability distribution of 
Qiviono in Fig.|| by using 3000 independent configurations at (3 = 2.4 on a 16 4 lattice. Several 
dotted lines correspond to partial contributions to the whole distribution, which are assigned 
to some integer value of Qsu(2)- We find several discrete peaks in Fig|| since each partial 
contribution is a Gaussian-type distribution with the half width slightly less than unity 
around discrete values (~ 0.7Qsu(2))- This implies that Qmoiio is classified by approximately 
discrete values. Namely, it seems that monopoles almost inherit the topological nature of 
the original gauge configurations. 
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IV. SUMMARY AND CONCLUSION 



Our work was motivated by the Ezawa-Iwazaki conjecture: the topological feature is 
preserved by the presence of QCD-monopoles once abelian dominance is postulated. To 
confirm this, we have performed an analytical study of how the monopole current is related 
to the topological charge within SU(2) lattice gauge theory. We consider the abelian analog 
of topological density, which is defined through the replacement of the SU(2) link variable 
by the U(l)-like link variable. Although this quantity is represented by the corresponding 
U(l) field strength in the naive continuum limit, it can carry the non-trivial contribution to 
the topological charge because of the presence of QCD-monopoles. If we assume that the 
abelian component of gauge fields plays an essential role in the description of the long range 
physics, QCD-monopoles are required to inherit the topological feature. 

Monte Carlo simulations have brought us encouraging results for the above conjecture. 
We have chosen the MA gauge as an applicable abelian gauge. We have measured both the 
ordinary topological charge Qsu(2) and the corresponding topological charge via monopoles, 
QMono- In spite of the fact that there is no correlation between Qsv(2) and Qtviono without 
abelian gauge fixing, we found a one-to-one correspondence between the two charges in the 
MA gauge. And also the relation: QMono ~ 0.7Qsu(2) can be read off from our resulting data. 
Furthermore, Qiviono is classified by approximately discrete values. Thus we have concluded 
that the topological nature is substantially inherited by QCD-monopoles. This conclusion is 
consistent with our previous study |33j , which shows that QCD-monopoles strongly correlate 
with the presence of fermionic zero-modes in the MA gauge. 
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FIGURE CAPTIONS 



FIG.l. The probability distribution of the amplitude of the off-diagonal compo- 
nents % on a 16 4 lattice at (3 — 2.4. 

FIG. 2. The topological density as a function of z and t in a two dimensional 
slice at the center of the instanton after 30 cooling sweeps (a), and after 
100 cooling sweeps (b) at (3 — 2.4 on a 16 4 lattice. 

FIG. 3. The abelian analog of topological density as a function of z and t in a two 
dimensional slice at the center of the instanton after 30 cooling sweeps 
(a), and after 100 cooling sweeps (b) at (3 — 2.4 on a 16 4 lattice. 

FIG. 4. Scatter plot of Qsu(2) vs. Qmoho (a) before the MA gauge fixing and (b) 
after the MA gauge fixing at 100 cooling sweeps. 

FIG. 5. Probability distribution of Qtviono i n the MA gauge fixing at 100 cooling 
sweeps by using 3000 independent configurations. 

TABLE 1. Slope of the linear correlation in the scatter plot of Qsu(2) vs. Qmoiio on 
a 16 4 lattice at (3 — 2.4, 2.45 and 2.5 after 30, 50 and 100 cooling sweeps. 
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FIG.l (Phys.Rev.D) Shoichi Sasaki et at. 




FIG. 2a (Phys.Rev.D) Shoichi Sasaki et al. 




30 cooling sweeps 



Fig. 2(a) 
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FIG. 2b (Phys.Rev.D) Shoichi Sasaki et al. 




100 cooling sweeps 

Fig. 2(b) 
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FIG. 3a (Phys.Rev.D) Shoichi Sasaki et al. 




30 cooling sweeps 

Fig. 3(a) 
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FIG. 3b (Phys.Rev.D) Shoichi Sasaki et al. 




100 cooling sweeps 



Fig. 3(b) 
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FIG. 4a (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. 4b (Phys.Rev.D) Shoichi Sasaki et al. 
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Fig. 4(b) 
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FIG. 5 (Phys.Rev.D) Shoichi Sasaki et at. 
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